Anomalous percolating properties of growing networks 
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We describe the anomalous phase transition of the emergence of the giant connected component in 
scale-free networks growing under mechanism of preferential linking. We obtain exact results for the 
size of the giant connected component and the distribution of vertices among connected components. 
We show that all the derivatives of the giant connected component size S over the rate 6 of the 
emergence of new edges are zero at the percolation threshold be, and S oc exp{— d(7)(& — bc)~^^^}, 
where the coefficient d is a function of the degree distribution exponent 7. In the entire phase 
without the giant component, these networks are in a "critical state": the probability 'P{k) that a 
vertex belongs to a connected component of a size A: is of a power-law form. At the phase transition 
point, T'{k) ~ l/(fclnfc)^. In the phase with the giant component, V{k) has an exponential cutoff at 
kc (X 1/S. In the simplest particular case, we present exact results for growing exponential networks. 
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I. INTRODUCTION 



From a physical point of view, networks may be equi- 
librium and non-equilibrium. For example, to the class 
of equilibrium networks belong classical random graphs 
with randomly distributed connections introduced by 
Erdos and Renyi and their generalizations |^|-^. 

Percolating 
studied lH,f 



roperties of equilibrium networks are well- 
-13|. The behavior of these networks near 
the threshold point, that is, near the point of the emer- 
gence of the giant connected component, is similar to 
percolation on an infinite-dimensional lattice. 

From the other hand, the most important real net- 
works (the Internet and the World Wide Web, for in- 
stance) have the growing total numbers of the vertices 
and, thus, are non-equihbrium Jl^-p^. The growth of 
networks, which is often a self-organization process, pro- 
duces a number of intriguing effects [|l5|,p0 21 1. 



Very recently, it was found that the percolation tran- 
sition in growing networks is of a quite different nature 
than in equilibrium ones [Q. (Note that the notions of 
percolation and a giant connected component are mean- 
ingful only in the large network limit.) In Ref. [^, for 
the growing network with an exponential degree distri- 
bution (degree is the number of connections of a vertex) , 
it was found numerically that this transition is of an in- 
finite order. All the derivatives of the size of the giant 
connected component (the percolating cluster) are zero 
at the percolation threshold. 

Here we propose a theory of the anomalous percolation 
transition in growing networks including the most inter- 



esting and important growing scale-free networks. Also, 
as a particular case, the exponential growing networks 
are considered. Thus, we present the complete exact de- 
scription of the percolation transition both for the scale- 
free and exponential growing networks. These cases have 
been turned to be similar to each other, and the phase 
transition is of an infinite order. Furthermore, we show 
that in the entire phase without the giant connected com- 
ponent such growing networks are in a "critical state". In 
this state, the probability 'P{k) that a randomly chosen 
vertex belongs to a connected component of the size k is 
of a power-law form. 



II. MODEL AND DEGREE DISTRIBUTION 

For constructing the scale-free network, we apply the 
mechanism of preferential attachments of new edges 



|l5| , p3| -|2q]. Here we use one of the simplest models pro- 
ducing power-law degree distributions 1 2^ : 

(i) At each increment of time, a new vertex is added to 
the network, so that the total number of vertices in the 
network is t. 

(ii) Simultaneously, b new undirected links are dis- 
tributed between vertices according to the following rule. 
The probability that a new edge connects a pair (/i, v) 
of vertices is proportional to (g^ -I- a)(g^ -I- a), where 
and Qi, are the degrees of these vertices, and a is some 
positive constant which plays the role of additional at- 
tractiveness of vertices for new edges b is also an 
arbitrary positive constant. Multiple edges are forbidden 
(in principle, they are non-essential for large networks). 
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The degree distribution P{q) of the network (degree is 
the total number of connections of a vertex) can be easily 
obtained using standard considerations (for example, see 



Refs. 



Pi'l) = 



33 



). It is of the form 



2bJ T{a) 



T{q + a) 



— ) 

2b I 



(1) 



where r( ) is the gamma-function. For large q, P{q) oc 
where 7 = 2 + a/ (26). It is convenient to introduce 
a new notation Q = ^ — 2 — a /(2b). In the limit — > cx), 
the 7 exponent approaches oo, and one can check that 
the degree distribution turns to be exponential. 



III. EVOLUTION OF CONNECTED 
COMPONENTS 

According to the above rules, a new vertex may have 
no connections, so disjoint vertices and components are 
certainly present in the network. We focus on the distri- 
bution of the sizes of connected components (clusters of 
mutually connected vertices) and on the size of the giant 
connected component. How do finite connected compo- 
nents grow with time? Here we present an elementary 
consideration of the evolution of connected components. 
For a rigorous derivation of our main equations, see Ap- 
pendix |A|. For the large network, it is almost impossible 
that both the ends of a new edge are being attached to 
the same finite connected component. Let us use the fol- 
lowing essential circumstance. In the large network,the 
finite size connected components are almost surely trees. 
Obviously, this is not the case for the giant connected 
component.) This fact is the key issue of percolation 
theory for networks |§,||J|] . One can check that the num- 
ber of edges in a finite connected component (tree) with 
k vertices is equal to fc — 1. Therefore, the total degree of 
this component equals 2(fc— 1), and the probability that a 
new edge is attached to the component is proportional to 
2k — 2 + ka. This value should be normalized. Taking into 
account that the total degree of the network is equal to 
2bt, we find that this probability is [{2+a)k-2]/[{2b+a)t]. 

The resulting equation for the number Nk{t) of con- 
nected components of size k at time t looks like 



Nk{t + l)^5{k-l) + 



1 - 26 



(2 + a)fc-2 



{2b + a)t 



Nk{t) 



(2 + a)j-2 (2 + a)(fc-j)-2 



(26 -f a)t 



{2b + a)t 



N,{t)Nu^,{t), (2) 



where 5{k— 1) is a convenient notation for the Kronecker 
symbol. The first term on the right-hand side of Eq. (|^) 
accounts for the emergence of new vertices. The second 
term describes the decrease of the number of connected 
components of the size k due to attaching of new edges 
to the vertices of these components. The last term on 
the right-hand side of Eq. (0) accounts for the fusion 



of connected components into larger ones of the size k. 
Here again we have used the largeness of the network to 
present the terms of the sum in the factorized form (see 
Ref. [^). The limit a ^ oo corresponds to the absence 
of any preference. In this case, new edges connect the 
pairs of randomly chosen vertices, and Eq. (2|) takes the 
form of the master equation derived in Ref. [ 22| . We em- 
phasize that Eq. (H) is nonlinear unlike master equations 
for degree distribution pqj29||. 



IV. MAIN EQUATIONS 

When t — > oo, the ratio Nk{t)/t approaches the sta- 
tionary value rifc. The equation for it is of the form 



1-1-26- 



26 



Uk = S{k- 1) -t- 



26 



2 



n-jUk-j, (3) 



Our main matter of interest is the probability V{k) that 
a randomly chosen vertex belongs to a connected com- 
ponent with k vertices, that is V{k) = kuk- Let us in- 
troduce its Z-transform: g{z) = X]fc°=i knkZ^ . Similarly, 
n{z) = X]fc°=i "^kZ^ ■ To approach our aim, that is the de- 
scription of the connected component statistics and the 
giant connected component, we must find giz). One sees 
that o(z) = zri'(z) and n(z) = dzg{z)/ z. Then, from 
Eq. (^) (see also Appendix ^), we obtain the equation 

, ^ 2 + a / , 2 , A 

g{z) - z + 26-— zg'iz) - ^r—g{z) x 
26 + a\ 2 + a J 



26- 



9{z) 



-n(z) 



(4) 



If the giant connected component is absent, that is, 
g{l) = 1, all the connected components are almost 
surely trees, so + o)k — 2]Nk = (26 + a)t, i.e., 

{2 + a)g{l) — 2n{l) = 2b + a, and we obtain the necessary 
condition: n{l) = 1 — 6 when 17(1) — 1. One can easily 
check that Eq. (|) satisfies this condition. 
Using the convenient combination 



hiz) ^ 1 - ^giz) 



dz 



26 + a 26 + a ./n z 

we rewrite Eq. (^ in the form 

d 



g{z) -z- bz^h'^iz) = 0. 
dz 



giz) , (5) 



(6) 



Therefore, 

n{z) 



dz 
z 



giz)=b[h^{z)-l] + z. 



(7) 
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Substituting Eqs. (^) and into Eq. (||), we finally 
obtain the equation 

(2 + a)zh{z)h'{z) - h^iz) + (1 + C,)h{z) - ((1 - z) = , 

(8) 

where C = 7-2 = a/(26). One sees that g(0) = n(0) = 0, 
so Eq. (||) must be supplied with the boundary condi- 
tion h{0) = 1. From Eqs. (H) and (0), a simple relation 
between g{z) and follows: 



obvious. Indeed, when C, = a/{2b) is fixed, from a ^ 
follows & — > 0. In turn, zero input flow of edges produces 
a set of disjoint vertices. 

For comparison, in Fig. 0, wc show the percolation 
threshold line for the equilibrium random graph with the 
same degree distribution (]l|) as our growing scale- free net- 
work [see the dashed line a(C) = (C — 3)/2]. This follows 
from the MoUoy-Reed criterion for the existence of the gi- 
ant connected component in equilibrium random graphs. 



2 + a 



l + lz+h^{z)~^h{z) 



(9) 



From Eq. (||), using Eq. (^, we can obtain the Z- 
transform g{z) of the basic distribution V{k) for finite 
connected components. When the giant connected com- 
ponent is absent, g(l) — 1 and n{l) = 1 — 6 as shown 
above. Then, from Eq. (||), it follows that h{l) = 0. As 
usually if the giant connected component exists, its 
size iS* can be related to .g(l): S = 1 — g{l)- After the 
introduction of new variable and function, y = z~^^'^'^^°^ 
and ip{y)/y = h{z)/{l + (), Eq. takes the canonical 
form of the Abel equation of the second kind 



C 



i + C 



(r 



y)- 



(10) 



The boundary condition for Eq. (|1C|) corresponding to 
the condition h{0) — 1 is^y) y/{l + as y oo. 

Setting ?/ = 1 in Eq. (Eo|), we see that two situations 
are possible. If (p{l) = 0, the giant connected component 
is absent. When (p{l) ^ 0, the giant connected compo- 
nent is present, and = 1. Then, one can obtain 
from Eqs. (||) or the expression for the size of the 
giant connected component 



^-1-5(1) = 



(1 + C)^ 



2 + a C 



^(1)[1-^(1)]. (11) 



Thus, our problem is reduced to the analysis of the solu- 
tions to Eq. (|0|). 



V. PHASE DIAGRAM 

The simplest problem we must solve is to indicate the 
region of the parameters b and a or, equivalently, C and 
a where the giant connected component is present. The 
direct analysis of Eq. ( p^ (see Appendix ^ yields the 
following picture (see Fig. |). When C> C = ^ + 2^2, 
the giant connected component is absent below the phase 
transition line 



(12) 



For ^ < C*, the trivial phase transition line is a = 0. The 
absence of the giant connected component at a = is 



VI. CRITICAL BEHAVIOR 

A. Size of the giant component 

In equilibrium networks, the size of the giant connected 
component linearly approaches zero at the phase transi- 
tion line. In growing networks we have a quite different 
situation. As shown in Appendix ^ the size of the giant 
connected component near the phase transition line takes 
the form 



C(2 + a) _ 1' 
(1 + C)2 4 



-1/2^ 



(13) 



The dependence of the factor D on ( is plotted in Fig. ||. 
One sees that, in the case of the network with an expo- 
nential degree distribution, i.e., when ^ — > oo, the factor 
D tends to a constant value 0.590 . . .. On the other side, 
D linearly approaches zero at C = C* ■ 

In the following, to simplify our expressions, we shall 
use the notation 



C(2 + a) _ 1 
(1 + C)2 4 



1/2 



(14) 



For the exponential network, i.e.. when ^ ^ oo, we have 
A= v/2|fo-l/8|.^ 

All the derivatives of the giant connected component 
size over the deviation from the critical line are zero. In 
particular, let us consider the deviation of b from the 
critical line be = &c(C) [the form of 6c (C) follows from Eq. 
(O)]. In this case, ^(C, b) is of the form 



SiC,b)^D{C)exp 



2V2 Vb~nrc 



(15) 



In the limit of the exponential network, —^ 00, where 
be = 1/8, we obtain 



^(C^oo,^) = 0.590... exp 



1 



2V2Vb~Fe 



(16) 



In Appendix we present a simple direct derivation of 
Eq. (1^). The factor in the index of the exponent is 
tt/ (2V2) = 1.111 .. . that is in agreement with the result 
of numerics in Rcf. |22| . 
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For small a and < C < C*j that is, near the other 
transition hne (a = 0, < C < C*) on the phase diagram, 
the size of the giant connected component behaves in the 
foUowing way 



S{C,a)^F{C)c 



(17) 



The factor F versus C is plotted in Fig. ^. All the deriva- 
tives of F over ( are zero at the point C = CI; Near this 
point, when a is small and = (* + e, Eq. (|l5| ) takes the 
form 

£"(£,0) = 0.631... aexp|-\/27r[a-4(3\/2-4)e]"i/2| _ 

(18) 

Equation (|l|) is valid when a-4(3\/2-4)e > 0. Thus we 
see that the phase transition of the emergence of the gi- 
ant connected component in growing networks is in sharp 
contrast to that in equilibrium nets. We shall disc uss t he 
nature of this anomalous phase transition in Sec. VII. 



We close this subsection by the exact result for the rel- 
ative total number £ of loops in the network (the ratio 
of the total number of loops in the network to its size 
t) obtained in the Appendix |A[ Note that all the loops 
are in the giant component since finite components are 
tree-like. The equation for £ looks like 



£=lj^[{2 + a)S + 2e]\ 



(19) 



Here we do not present its trivial solution. Near the 
phase transition, where the giant component is small, 
£-[(2 + a)5]V[2(l + C)]. 



when the parameter A is not small [see the definition 
©I, 



r{k) ^ fc-4A/(l-2A) 



(21) 



In the same phase, near the phase transition line (p^), 
V{k) has a power-law tail 



32a 



X^[{2 + a)k] 



-2-4A 



(22) 



for Ink ^ 1/A, and coincides with the threshold distri- 
bution (|20|) in the region 1 < Infc < 1/A. 

This power-law form is in striking contrast to the ex- 
ponentially decreasing V{k) both above and below the 
percolation threshold for standard percolation in- 
cluding percolation on equilibrium networks. Equation 
( pO| ) indicates that the growing network is in a "critical 
state" in the entire phase without the giant connected 
component. Note that this is valid both for the scale-free 
and exponential growing networks. 

In the phase with the giant connected component, 
V{k) has an exponential tail. Near the phase transi- 
tion, for the large values of connected component sizes 
k ^ [a/ (2 -I- a)]S~'^, we obtain the following behavior: 



■P(fc) [27r(2ae)'^(2 + a)5-'^fc] 



-3W-1/2 



exp[ 



2ae 



■Sk]. (23) 



Here, a — a{() is taken for the transition line (^^, and 
S — S{a, is the size of the giant component given 
by Eq. (p^; e = 2.73.... For smaller fc, namely 
1 ^ In/c <C 1/A, 'P{k) takes the form of the threshold 
distribution (EO). For Infc ~ 1/A, a crossover regime is 
present. 



B. Distribution of vertices among connected 
components 

The distribution of vertices among connected compo- 
nents is one of the key issues in percolation theory. In 
Appendix we calculate the probability V{k) that a 
randomly chosen vertex belongs to a connected compo- 
nent of the size k. 

At the point of the emergence of the giant connected 
component, we obtain 



VII. INTERPRETATION 



Vik) 



2a 



1 



2 + a fc2 In^ /fc 



(20) 



Here, a — a(C) is given by the relationjlj) for the critical 
line. The factor 2a/(2 -I- a) in E n. ( [2C ) approaches zero 
at the point C = C*- Equation ( |20| ) must be compared 
with the corresponding result for percolation on the equi- 
librium networks and infinite-dimensional lattices, where 
V{k) cx k-^/^ at the percolation threshold |3[|,|,||l. 

Furthermore, we find that in the entire phase without 
the giant connected component, V{k) is of a power-law 
form. In this phase, far from the phase transition, i.e.. 



All the results of Sees. ^ and VI were obtained by 
the explicit formal analysis of the master equations for 
growing networks. Let us discuss the physical nature of 
the behavior observed above. 

To begin with, let us recall that in networks growing 
under mechanism of preferential attachment of edges to 
vertices degree distributions of vertices are of a power- 
law (fractal) form. This means that such networks self- 
organize into scale-free structures with the power-law de- 
gree distributions while growing. In fact, they are in a 
critical state in a wide range of the values of the network 
parameters. The growth under mechanism of the prefer- 
ential attachment produces power-law distributions. In 
principle, this phenomenon can be called self-organized 
criticality. 

In the present paper, we are interested not in the statis- 
tics of vertices but in the statistics of connected compo- 
nents, that is, in the distributions of the number of con- 
nections of distinct connected components and the dis- 
tributions of the number of vertices in them. New links 
are being attached to large connected components with 
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higher probabihty, so that large connected components 
have a better chance to merge and grow. This produces 
the preferential growth of large connected components 
even in networks where new edges are attached to ran- 
domly chosen vertices, that is, in networks without pref- 
erential attachment of edges to vertices. Such mechanism 
of the effective preferential attachment of new vertices to 
large connected components naturally produces power- 
law distributions of the sizes of connected components 
and power-law probabilities 'P{k). This "self-organized 
critical state" is realized in the growing networks only if 
the giant component is absent. 

As soon as the giant connected component emerges, 
the situation changes radically. A new channel of the evo- 
lution of the connected components is coming into play, 
and, with high probability, large connected components 
do not grow up to even larger ones but join to the gi- 
ant component. Therefore, there are few large connected 
components if the giant component is present, and then 
V{k) is exponential. 

Thus, in the growing networks, two phases are in con- 
tact at the point of the emergence of the giant connected 
component — the critical phase without the giant com- 
ponent and the normal phase with the giant component. 
This contact provides the above observed effects. There 
exists another example of a contact of a "critical phase" 
(or of a line of critical points) with a normal phase, 
namely, the Berezinskii-Kosterlitz-Thouless phase tran- 
sition |35|j3q| . Interestingly, functional dependences in 
both these cases have similar functional forms. This in- 
dicates that equations describing such phase transitions 
have similar analytical properties. 

Connections in non-equilibrium networks are very in- 
homogeneously distributed between vertices. We mean 
that many edges are captured by old vertices, and 
few edges are attached to more young vertices — "the 
rich gets richer". Nevertheless, this statistically in- 
homogeneous distribution of connections in growing net- 
works is not the direct origin of the observed behavior. 
Both this in-homogeneity and the critical (or, one can 
say, power-law, or fractal, or scale-free) distributions of 
connected components in the absence of the giant com- 
ponent have the same first cause — the specific process 
of the network growth. 



VIII. CONCLUSIONS 

In summary, we have presented the theory of percola- 
tion in evolving systems — growing networks. We have 
demonstrated that the interplay of a self-organization 
process and percolation produces a number of intriguing 
effects in such objects. We have obtained exact results 
for the size of the giant component and the distribution of 
vertices over connected components. An explicit descrip- 
tion for the anomalous phase transition of the emergence 
of the giant component in these growing networks has 



been proposed. We hope that our results are of a general 
nature and can be applied to various growing systems. 
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APPENDIX A: RIGOROUS DESCRIPTION OF 
THE EVOLUTION OF CONNECTED 
COMPONENTS 



In Sees. Ill and [V , we have derived our main equation 
(?? using a simple but rather heuristic approach. Here 
we present a strict derivation of the equations describing 
the growth of the network. 

Let N{t) be the number of vertices in the growing net- 
work at time t. We assume that, with the probability 
p(t)dt, a new vertex is added to the network during a 
small time interval dt, i.e., M{t + dt) = Af{t) + 1. The 
degree of a new vertex is supposed to be zero. 

The total degree of the network is 



Q{t) = 2C{t) 



M(t) 



(Al) 



i=l 



where C{t) is the total number of edges in the network. 
We assume that with probability b{t), a new edge emerges 
between vertices. According to the rule of preferential 
linking that we use in the present paper, the probability 
that this edge connects vertices i and j is equal to 



[qi{t) + a][qj{t) + a] 
[Q{t) + aN{t)Y 



(A2) 



for each pair of vertices. 

If h{t) = const, one can set p{t) — const — 1, so, in this 
case, we have two parameters, namely b and a — addi- 
tional attractiveness, that determine the growth and the 
structure of the network. Let us introduce a new object 
which we call the connectivity matrix of the network: 

{1, if i and j belong to the same 
connected component (A3) 
0, otherwise 

(this should not be mixed with the adjacency matrix). 

Then ki It) — Y^^l^i Stj (t) is the size of a connected com- 
ponent containing i-th vertex, that is, the total number 
of vertices in this component. 
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The closed equation can be written for the joint dis- 
tribution function of the total number of vertices in the 
network, M, the total degree, Q = 2£, and the connected 
component sizes, {ki}. Let us define the Z-transform of 
this joint distribution function as 

g{U, Q; x,t)^^ (s[m ^] S[Q{t) - Q] ^ x'^'A . 



Recall that (5( ) is the Kronecker symbol. Then, 
g{Af,Q;x,t + dt) = {1 - dt - bdt)g{N', Q;x,t) 



(A4) 



dt 



J2 (Sl^it) + 1 - AA] 5[Qit) - Q]x'=' 



(«) 



bdt 

w 



{5[^^{t)-^^]5[Q{t)^2-Q]x 



[q,{t) + a][qi{t) + a] 



AT 



_ ^fcj(*)+fel(t)l 



.m— 1 



m — 1 



(A5) 



In the large Af, Q limit, one can substitute the factor 
[1 - Suit)] in Eq. (^) for 1, therefore 



.g(7V, Q;x,t + dt) ^ [1 - dt -bdt)g{N, Q; x, t) + 

6[Af{t)-Af+\\5[Q{t)-Q]x''^^''> 



N N- 1 



^{5[M{t)-N+l]m{t)~Q])x- 
bdtg{Af,Q-2;x,t) - 
2b dt 



JV{Q-2 + aM) 



{d[Af{t)-Af]S[Q{t)-Q + 2]x 



Y,[q,{t)+a]k,{t)x'^^^') 



2b dt 



N{Q - 2 + oNf 



{5[N{t)-N]5[Q{t)- Q + 2]-K 



(A6) 



N(t) M(t) N(t) 

i=i i=i j=i 

E x'^^W E - 2 E - 1]^'^-^*^ ■ (A7) 

j=l 1=1 j=l 

Here we have used the fact that, in the tree ansatz, loops 

are absent, so Yl^=i = 2[fcj(t) — 1]. Analo- 
gously, 



A/-(t) 



A/-(t) 



4=1 



rfej (t) 



(A8) 



One can see that a general relation 



E<z.w/[fc4W] = 2E 

i=i j=i 



h{t) 



mm m 



holds for arbitrary graphs with tree-like finite compo- 
nents. Here /( ) is an arbitrary function. In particu- 
lar, if f{k) = 1, then 2Q{t) = 2[Ar{t) - Afcom{t)], so 
J^{t) — Qit) +Ncom{t), where Mcom{t) is the number of 
components in the network at time t. Using these rela- 
tions, we obtain 

g{Af,Q;x,t + dt) = {1 - dt - b dt)g{Af, Q; x,t) + 



l-^)g{N-l,Q; x, t) + ^g{Af - 1, Q; 1, t)x - 



bdtg{M,Q-2;x,t) 
2b dt 



Af{Q-2 + aAf) 
Y}{'^ + a)h{t)-2]i 



{5[n{t)~n]6[Q{t)-Q + 2]x 



Mt) 



2b dt 



n{Q-2 + an)'' 



{5[Af{t)~U]5[Q{t)- Q + 2]-K 



E[(2 + a)A:,W-2]x'=-«||E 

3 = 1 



2 + a- 



h{t) 



(AlO) 

In particular, when a; = 1, one obtains the joint prob- 
ability that the network contains A/" vertices and and 
C = Q/2 edges at time i, 

g{U, Q; 1, t) = {S[U{t) - Af] 5[Q(i) - Q]) = n(A^, Q; t) . 

(All) 



The sums in Eq. (A6) can be easily calculated: 



From Eq. (AlO), we obtain the equation for n(A/', Q;t): 



6 



dUjU, Q; t) 
dt 



= Tl{U -l,Q;t)-n{U, Q;t) 



(A12) 



Choosing the initial condition Il{JV,Q;to ) = d{JV — 
J^o)S{Q - Qo), we find the solution of Eq. ([AT2|), 



n(AA,Q;t): 

^^f-^ro 
(AA-ATo)! 



^l + (-l)g-Qo (bt)(g-g°)/2 
2 [{Q-Qo)/2]\ 



(A13) 



If we assu me t hat Qq — 2btQ and A/q = <o for tQ ^ oo, 
then Eq. ( |A1^ ) yields n{Af, Q:t) 6{M - t)5{Q - 2bt). 
This shows that the total numbers of vertices and edges, 
M and Q, are, in fact, rigidly determined in the large 
network limit. Finally, using the decomposition 

(^fe-a^fe.) ^ (A14) 

which can be justified in the limit of t — > (X), from Eq. 
( |A1C| ), we obtain 



,d9_ 
'dt 



9 



2b 



2^ N dg 



-^[(2 + a).£ 



(26 + a) 



25]- 

2g][{2 + a)g - 272] = X , (A15) 



wher e n{x ,t) = Jq dyg(y,t)/y. In the stationary case, 
Eq. ( [aTsI ) yields Eq. (|) of Sec. |v. 

If the giant component is absent, = 
limj;^i_jV'-»oo 5(2;, A/"; = 1. In addition, if all the fi- 
nite connected components in the network are tree-like, 
we have 



n(l,t) = 



1 



E 



1 



,m hit)/ 

( M{t)-Q{t)/2 \ 

\ m) I 



zom (^) 



(AI6) 



that is, the "tree condition" . The following equation can 
be written for n{x,t): 



dn 



2b 



2b 



■[i2 + a)g-2n] 



{2b + a)' 



:[{2 + a)g - 2nf ^x, 



(A17) 



Equation ( |A17| ) can be used for the determination of 
the number of loops in the giant connected component 
which coincides with the total number A4 of loops in the 
network since all the finite connected components are 
tree-like. Ai indicates the extent of the deviation of the 
giant component structure from a tree. One sees that 
A^ = £ -I- 1 — Af. £ = bt,Af = t, hence the total number 
of edges in the giant component is equal to 



bt 



J2ik - l)Nk{t) = t[b - Y^ik - l)n,(t)] 



t[b-g{l) + nil)] 



(A18) 



Subtracting from this expression the total number of ver- 
tices in the giant component, that is, tS — t[l ~ g{l)], we 
obtain A^ = t[n(l) — (1 — 6)]. It is convenient to introduce 
the relative number of loops for i — > 00: 



^^:^ = n(l)-(l-6). 



(A19) 



Substituting n(l) from Eq. (A19) into Eq. (A17) taken 
at the point x ~ 1, we obtain the exact equation for £: 



2b + a 



[{2 + a)S + 2£f 



(A20) 



This equation also follows from Eq. (^). Near the phase 
transition, the giant component is small, so that 



^ (2 + a)^ 
' 2b + a 



(A21) 



APPENDIX B: ANALYSIS OF EQ. ^ 

Note, that Eq. ( p^ ) remains unchanged, if C ^ 
However, the condition that (p/y ^ 1/ (1 -I- C) at y ^ 00, 
is not invariant under this change, and for C < 1 we 
have to choose another solution then for C > 1- Here 
we analyse the solutions of Eq. (10) at y — > cx) and at 
?/ — !■ 1. The comparison of both the regions will allow 
us to obtain the phase diagram and to find the essential 
features of the probability P(fc) that a vertex belongs to 
a connected component of the size k. 

If y 00, the term y"^^'^ on the right-hand side of 
Eq. ( |lO| ) may be neglected, and we obtain the equation 



tptp-tp 



c 



(1 + 



(Bl) 



This equation has two solutions, which are linear in y: 
= + and (f = Cy/(1 + C)- It is the first 

one, which must be chosen, because it corresponds to 
g (0) — 0. Let us denote the physical solution of Eq. 
( pO| ) as (pi, and as (p2, — the other one, which has the 
asymptotic form ^y/ (1 -I- (^) at y — s- 00. If C > 1, then 
ifi (y) < ip2 (y) for y > 1. This follows from the unique- 
ness property of a solution of Eq. ([lo| ) . On the contrary, 
if C < 1, the physical solution is a higher one, (fi > (^2- 
At y 1, after the linearization of the right-hand side 



of Eq. (10) with respect to y — 1 we obtain: 



ipip-ip 



-+ P] (y-1) : 13 = 
4 (1 + C) 



^(2 + «)-7 



(B2) 
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After the substitution ip = {1 — y) , Eq. (B2) takes the 
form 



and can be easily solved: 
ln[C(2/-l)] = 



■0 dip 



Here C is the integration constant. The sign of /3 de- 
termines a full picture of the set of the solutions of Eq. 

(H- 

At first, l et u s consider the case /3 = — < 0. In this 
case, Eq. (B2) has three families of solutions, real for 
y > 1. These families can be written in the implicit form 
as: 



C(y-l) 



1-y 2 

Ciy-l) = 



2 1-y 



2 1-y 



_1 _1 

4A 2 



1 

4A 2 



1 

4A 2 



^ 1-A 



1-y 2 



^ i-A 



1-y 2 



+ A 



1-y 



_i 1 

' 4A 2 



' 4A 2 



1 

4A 2 



(B3) 



(B4) 



(B5) 



In the family (gaj), 1/2 + A < ip{l - y) < +oo, in (gj), 
1/2- A < 'p\l-y) < 1/2 + A, i n (|B5|) , -oo < fil-y) < 
1/2 — A. Only the families (|B3|) and (|B5[) should be 

we have 



A. Only the families (p3|) and (|B5|) 
taken into account, because, for the family (| 
g" ( 1) < 0. Here a physical solution is r ealized in the fam- 
ily dBl) if C < 1, and in the family (|B5h when C > 1. The 
distinctive feature of the solutions (B3) is that they have 
nonzero value asy^l: (p{l) — 1 JC. This proves that 
when C < 1: the giant com pon ent is always present. For 
the solutions of the family (B5) we have (1) = 0, which 
means the absence of the giant component for C > 1 and 
/3 < 0. When y ^ 1, we obtain from Eq. (B5): 



\-^ + C'{y 



1) 



4A/(1-2A) 



(B6) 



where C is some constant of the order of unity if A ~ 1 . 

If /3 = A^ > 0, the solution may be conveniently ex- 
pressed as: 



C(y-1) = 



1-y 



+ A2 



-1/2 



This form of presentation wa s chose n to ensure a smooth 
crossover between Eqs. ( |B5| ) and (B7) at small A. All 
solutions of this set take nonzero values as y — + 1: 



^(1) 



1 

C 



exp 



(B8) 



This means, that for C > 1, the value /3 = corresponds 
to the point of the emergence of the giant connected com- 
ponent in th e ne twork. Taking into account the definition 
of /3 in Eq. (B2), we arrive at the expression jl^ ) for the 
critical line in the a) plane. 

Let as consider now the critical region, |/3| ^ 1. If 
/3 = 0, the solutions of Eq. (B2) are given by 



y-1 



1 



1 



W[-Ciy-1)] 



(B9) 



where, by definition, the Lambert function W (z) is a 
proper solution of the equation W exp W = z. If —e~^ < 
z < Q,W has two real branches: the one for which W ^ 
as z —0, and the other for which W — > — oo in the 
same limit, z = —e^^ is the branching point for both 
these branches. When the integration constant C is pos- 
itive, one must choose that real negative branch of W (z), 
which tends to — oo as y ^ 1 -f 0. This ensures a smooth 
crossover between Eqs. (B5) and (B7) as A ^ 0. 

The behavior of the distribution near the critical line 
may be treated analytically, if w e know the integration 
constant C — C {(,) in Eq. (B9), that may be obtained 
by numerical integration of Eq. ( [lO| ) at the critical line, 
given by Eq. (12). Indeed, as 1 3> y — 1 ^ exp (— 7r/2A), 
the argument of arctan in Eq. (B7) becomes nega- 
tive and large. Then, using the asymptotic expression: 
arctan z w — 7r/2 + 1/z, one can see that the solution 



(B7) turns into the solution (B9) with the same integra- 
tion constant C. Hence, substituting Eq. (B8) into the 
expression for the size of the giant component (pi]), ac- 
counting for the relation ( |l2|) and recalling that A = \/]3, 
where /3 is defined in Eq. (p2[) , one arrives at the expres- 
sion ( |l3|) for the giant component size. The resulting 
factor D{() in this expression is equal to 



^(0 = 



2a (C) ^ - 6C + 1 

c(c) 2CC(C) 



(BIO) 



Now let us consider the distribution function for con- 
nected components V (fc) at the threshold and near it. 
This is the inverse Z-transform of g (z): 



V{k)=<p^^g{z)z 



(Bll) 



where the integration is performed along the contour 
around z = 0, lying inside the unit circle. After inte- 
gration by parts, accounting for Eq. (^, this expression 
takes the form 



exp ■ 



1 

'2A 



arctan 



A Vl-y 



(B7) 



V (fc) = S{k-l) + ^f (z) z~^-' . (B12) 
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Introducing the integration variable y = x 
(p{y), ^Iv = h/ (1 + C) we obtain: 



-l/(2+a) 



V{k) = 5{k~l)+a{2 + a) 



(1 + C)^ 
2C 



(2+a)fc 



where c is some integration contour, lying to the right of 
y — 1 point. At /c 3> 1, the position ?/c and character of 
singularity with the highest value of jt/d determine the 
value of this integral. Close to the transition line we have 
either the singularity at j/c = 1, if /9 < 0, or at ?/c = 1 — 
e <C 1, if /? > 0. Hence, when A = y/P :s> 1, and k is 
large enough, the vicinity of y = 1 yields the main con- 
tribution to the above integral. Then one can extend the 
integration contour to ±ioo. Changing the integration 
variable: y = 1 + s, and assuming that s is small, we 
finally obtain the expression for the large k part of the 
connected component distribution. 



r{k) w a(2 + a) 



(1 + Cf 
2C 



ds 
2ni 



<^2(l + s)e(2+-)'=^ 



(B13) 



At the threshold we have: 

„2 



In (Cs) 



where we have taken into account, that the appropri- 
ate branch of Lambert function has the asymptotic form 
W (— z) w \nz as \z\ <C 1. Then, deforming the in- 
tegration contour to the one along the shores of the 
s = (— oo, 0) cut, and calculating the jump across the 
cut to the leading order in 1/lnfc, we obtain: 



Vik) 



1 



2 + ak'^ln^k 



(B15) 



When /? < 0, that is, in the phase without the giant 
connected component, and the strong inequality A <C 1 
is not valid, we have from Eq. (B6): 



^(1 



A +Ais^ 



4A 



2A' 
(B16) 



with some constant coefficient ~ 1. Only the second 
term in the square brackets is singul ar an d yields no nzero 
contribution. Substitution of Eq. (B16) into Eq. ( B13| ) 
yields: 



r{k) « At 



-2-A 



(B17) 



which is valid if In [(2 + a) k] l/X; A is a constant. At 
smaller k one should use expression for ip(l + s), valid 
at larger s. If Ins ^ 1/A, from Eq. ( p5| ) the equation 
ioi (p{l + s) follows: 



Cs = 



exp ■ 



(B18) 



whose solution is given by the fun ction ip{y), taken pre- 
cisely at the threshold, Eq. (^9|). Therefore, as 1 ^ 
In [(2 -I- a) fc] ^ 1/A, the distribution function assumes 
the threshold form, Eq. (B15). 

Above the threshold (i.e., in the phase with the gi- 
ant connected componen t), > 0, the argument of the 
arctan function in Eq. (B7) is positive and large, and 
the formula arctan z«7r/2— 1/z may be used. Thus, we 
obtain the form of the solution: 




1 



W[Csexp(^)] 



(B19) 



This expression is valid if ln[sexp(^)] ^ 1/A. The 
Lambert function W {z) has a square root type sin- 
gularity at z = — e —1, which ensures the exponen- 
tial type behavi or of V (fc) at th e lar gest k. Let us 
substitute Eq. ( p319| ) into Eq. (B13), retaining only 
the relevant term. Changing the integration variable, 
s = C~^zexp (— ^ + z), we find the expression for the 
distribution: 



(B14) Vik) « C'-'a (2 + ay /c exp - 



2A 



dz 
27ri' 



z (z + 1) exp [C^^ (2 + a) e^'^/^^fcze^ 3zl . (B20) 



If {2 + a) k ^ exp(^), this integral can be calculated 
in the saddle point approximation. Notice, however, that 
the integrand in Eq. (B20) becomes zero at the saddle 
point Zc = — 1. To avoid this difficulty, one can perform 
integration by parts, which gives: 



A 



dz 
27ri 



(2z -t- 1) exp [C^i (2 -t- a) e^'^/^^fcze^ 2zj . (B21) 



Then, the saddle point approximation yields: 
P (/fc) « a \2ttC^ (2 + af k^ 



exp 



(B22) 



At smaller 
less, it is 
i.e., when (2 
|ln [(2 -I- a) A: exp (- 
W-function in Eq, 



values of k 
still possible to 
a) k exp ( 




but when, neverthe- 
use Eq. 

< 1 



2\/ 



( IB19D , 
but 



^ 1/A, the argument of the 
19() becomes large, and we can re- 
place the Lambert function with logarithm. In the same 
way as obtaining the threshold distribution (B15), we 
get the form of the distribution: 
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V{k) 



2a 



(2 + a)fcexp(-^)] . (B23) 



At even smaller k, the expression for (1 + s) at larger 
|s| is necessary. It may be found from Eq. (B7), if we 



assume that the argument of the arctan function is small 
and replace arctan z with z. As a result we obtain 



1 - In ( 2\Cse''''^^ 



(B24) 



which is valid if |ln [Xse"/^^^ | <C 1/A. In this region the 
distribution function is of the form 



V (k) w 



(B25) 



Finally, for k ^ 1, \nk ^ 1/A, equation (B7 ) fo r 
tp{l + s) at s ~ l/[(2 + a)A:] assumes the form (B18), 
and as a result we obtain the distribution function in its 
threshold form, Eq. (Bl^). 



APPENDIX C: ANOTHER WAY TO GET S{B) 
FOR THE EXPONENTIAL NETWORK 

Here we show how our result for S{b) can be obtained 
directly for the network growing without preferential at- 
tachment of edges, i.e., in the limit a — s- cx). In this par- 
ticular case, from Eq. (^), we obtain the master equation 
for the probability V{k): 

dVik) 



dt 



fc-1 



5{k -\) + bkY^ V{i)V{k - j) - 2bkV{k) . (CI) 

i=i 

This is a basic equation for the evolution of the c onn ected 
components. From the long-time limit of Eq. (CI), the 
equation for g[z) follows (see Ref. |22[): 



If . 1 z- g{z) 



(C2) 



The boundary condition for it is g(0) = 0, so g'(0) = 
1/(1 + 2&). 

The threshold solution g{z,b = 1/8) approaches 1 at 
z = I, and g'{l,b = 1/8) = 2. For b > 1/8, the giant 
connected component is present, so that, at z = 1, the 
corresponding solution of Eq. ( |C^ ) is less than 1, and 
g'{l) = l/{2b). When 6 < 1/8, i.e., in the phase with- 
out the giant connected component, the physical solution 
g[z = 1) = 1, and g'{l) = (1 - Vl - 86)/46. Note that 
g{z) approaches the point z = 1 in a non-trivial way. In- 
deed, the values of g'{z) are essentially smaller than g'{l) 
even very close to z = 1 (see below). 

Near z = 1, Eq. (C2|) can be written in the form 



(C3) 



where ^ = 1 — z and u = \ — g. Its solution for b = 1/8, 
that is, the threshold solution, is 

^.(e,5-l/8) = 2e[l-/(0], (C4) 
where /(^) is the solution of the transcendental equation 

1 



Hum 



= In c . 



(C5) 



Here, the constant c — 0.295 ... is obtained by the numer- 
ical sewing together with the solution g{z) of Eq. (C2) 
passing through at z = 0. 

For b > 1/8, i.e., when the giant connected compo- 
nent is present, the solution of Eq. (C3) is given by the 
following transcendental equation 



1 



V8b^ 



arctan 



4b[u{0/^] - 1 
V8b~T 



(C6) 



One can check that du{0) /d^ — 1/ {2b) . The constant C is 
fixed by the value of the solution at ^ = 0, i.e., u(0) = S, 
where S is the size of the giant connected component: 



C 



tt/2 



\n\/2b-\nS. 



(C7) 



When b tends to 1/8 from above, for ^ 3> 5, we expand 
arctan [accounting that u(^) < 2^ in this region] and set 
& to 1/8 in the logarithms: 

^/^-mi-in^^ 

1 2 



V86 
7r/2 

V8fe^ " 1 - u{OI{2£) 



1 



InK - u{0/2] 



(C8) 



We sew together this solution for 6 ^ 1/8 and the thresh- 
old one (|C^) . One can see that this is possible substitut- 
ing Eq. (ra) into Eq. (|C8|): 



Accounting for Eq. (|C5|), we finally obtain 



S — 2c exp 



1 



2V2 ^fe- 1/8 



(C9) 



(CIO) 



where the coefficient 2c — 0.590 . . .. 

One should note that the accurate sewing procedure 
has been necessary only for the deter mination of the co- 
efficient of the exponent in Eq. ( plOj ). Indeed, the index 
of the exponent can be easily obtained without consider- 
ation of the last two terms on the right-hand side of Eq. 
(S for g{x). 
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FIG. 1. Phase diagram of the growing scale-free network. 
Here, a is additional attractiveness of vertices for new edges, 
and = 7 — 2 = a /{2b), where b is the value of the input flow 
of new edges. The solid line [see Eq. (|l^)] indicates the phase 
transition of the emergence of the giant connected component. 
When < 3 + 2\/2, the giant connected component is absent 
(5 = 0) only on the a = Q line. The dashed line shows the 
percolation threshold line for the equilibrium random graph 
with the same degree distribution (^ as the growing scale-free 
network under consideration. Inset: the same phase diagram 
but plotted on axes (C,^)- 
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